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£ry It appears that the majority of rapidly-, well-localized gamma-ray bursts with 

undetected, or dark, optical afterglows, or 'dark bursts' for short, occur in clouds of 
size R ^ 10L 4 g pc and mass M^3x 10 5 L4g Mq, where L is the isotropic-equivalent 
peak luminosity of the optical flash. We show that clouds of this size and mass cannot 
be modeled as a gas that is bound by pressure equilibrium with a warm or hot phase 
of the interstellar medium (i.e., a diffuse cloud): Such a cloud would be unstable to 
r^ gravitational collapse, resulting in the collapse and fragmentation of the cloud until a 

burst of star formation re-establishes pressure equilibrium within the fragments, and 

P the fragments are bound by self-gravity (i.e., a molecular cloud). Consequently, dark 

bursts probably occur in molecular clouds, in which case dark bursts are probably 
a byproduct of this burst of star formation if the molecular cloud formed recently, 
and/or the result of lingering or latter generation star formation if the molecular cloud 
formed some time ago. We then show that if bursts occur in Galactic-like molecular 
clouds, the column densities of which might be universal, the number of dark bursts 

i_> . can be comparable to the number of bursts with detected optical afterglows: This is 

what is observed, which suggests that the bursts with detected optical afterglows might 

£> \ also occur in molecular clouds. We confirm this by modeling and constraining the 

distribution of column densities, measured from absorption of the X-ray afterglow, of 



J2 



X 



the bursts with detected optical afterglows: We find that this distribution is consistent 
with the expectation for bursts that occur in molecular clouds, and is not consistent 
with the expectation for bursts that occur in diffuse clouds. Consequently, we find that 
all but perhaps a few bursts, dark or otherwise, probably occur in molecular clouds. 
Finally, we show that the limited information that is available on the column densities 
of the dark bursts is not consistent with the idea that the dark bursts occur in the 
nuclear regions of ultraluminous infrared/submillimeter-bright galaxies, from which we 
draw conclusions about the nature of star formation in the universe. 
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1. Introduction 

In Reichart & Yost (2001), we show that about 60% of all bursts, and about 80% of rapidly-, 
well-localized gamma-ray bursts with undetected, or dark, optical afterglows, or 'dark bursts' for 
short, have afterglows that are fainter than R = 24 mag 18 hours after the burst. Furthermore, we 
show that, with the exception of perhaps a few bursts, the dark bursts are most likely the result of 
circumburst^] extinction, and the density of the circumburst medium probably spans many orders 
of magnitude, from densities possibly as low as densities that are typical of the Galactic disk to 
densities probably as high as densities that are typical of dense clouds. 

In Reichart (2001b; see also Waxman &; Draine 2000; Fruchter, Krolik & Rhoads 2001), we 
show that sublimation of dust by the optical flash, and fragmentation of dust by the burst and 
afterglow, prevent clouds that are smaller than R ~ 10L 4 g pc, where L = 10 49 L4g erg s _1 is the 1 
- 7.5 eV isotropic-equivalent peak luminosity of the optical flash, from significantly extinguishing 
their afterglows. Given the finding of Reichart & Yost (2001) that circumburst extinction appears 
to be responsible for most of the dark bursts, we find in Reichart (2001b) that most of the dark 
bursts occur in clouds of size R ^ 10L 4 g pc and column density Nu ^5x 10 21 cm -2 , and hence 
mass M^3x 10 5 L 49 M . 

In this paper, we consider the equilibrium properties of clouds, and show that clouds of these 
sizes and masses, and hence the circumburst clouds of dark bursts, cannot be diffuse, but rather 
are probably molecular, in which case they are active regions of star formation (§2). In §3, we show 
that Galactic-like molecular clouds have sufficiently high column densities to make bursts dark at 
optical wavelengths, and that sublimation of dust by the optical flash and fragmentation of dust by 
the burst and afterglow can leave a comparable number of bursts - bursts that are sufficiently near 
the earth-facing side of the cloud and/or in sufficiently small clouds - relatively unextinguished at 
optical wavelengths, which is in agreement with what is observed. In §4, we model and constrain 
the column density distribution of the bursts with detected optical afterglows, and show that 
these bursts are also consistent with a molecular cloud origin, and are not consistent with a diffuse 
cloud origin. Consequently, we find that all but perhaps a few bursts, dark or otherwise, probably 
occur in molecular clouds. 

In §5, we discuss implications for the isotropic-equivalent peak luminosity of the optical flash, 
the density of the circumburst medium, circumburst extinction studies, and the nature of star 



By 'circumburst', we mean within the circumburst cloud. 



formation in the universe: We find that for most bursts, the isotropic-equivalent peak luminosity 
of the optical flash is probably less than L49 ~ 0.1 (§5.1), that the density of the circumburst 
medium is possibly higher than what efforts to model X-ray through radio afterglow data have 
found to date (§5.2), that bursts might serve as probes of dust under the extreme conditions of 
sublimation and fragmentation, and as a function of redshift (§5.3), and that our finding that 
most, if not all, bursts occur in molecular clouds strengthens the pre-existing idea that most of the 
star formation in the universe occurs in molecular clouds (§5.4). We summarize our conclusions in 
§6. 



2. Evidence for a Molecular Cloud Origin for Dark Bursts: Cloud Equilibrium, 
Gravitational Collapse, and Fragmentation-Driven Bursts of Star Formation 

We first consider the equilibrium properties of diffuse clouds: i.e., gas that is bound by 
pressure equilibrium with a warm or hot phase of the interstellar medium. Specifically, Spitzer 
(1978) considers the equilibrium properties of a uniformly magnetized, isothermal, non-rotating 
diffuse cloud, and introduces factors c\ = 0.53 and c 2 = 0.60 that correct for the fact that in 
equilibrium the cloud will not be of uniform hydrogen density tih, but will be centrally condensed, 
and will not be spherical with radius R, but will be moderately flattened along the direction of 
the magnetic field. Spitzer (1978) shows that clouds with R > R m are unstable to gravitational 

collapse, where 

GM 2 = 2hy m 
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M = AirpR^/S is the mass of the cloud, p = fj>H m H n H, fJ-H = 1-87 is the mean molecular weight 
per hydrogen atom, 
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and B is the strength of the magnetic field within the cloud. We plot R m {riH) for T = 10 and 100 
K and B = 3 pG in Figure 1, and for T = 30 K and B = 1 and 10 pG in Figure 2: Most Galactic 
clouds have 10 < T < 100 K, and B = 3 pG is typical of the Galactic interstellar medium. 

Also in Figures 1 and 2, we plot the curves of constant post-sublimation/fragmentation 1 - 
7.5 eV optical depth r from Figures 3 and 4 of Reichart (2001b): The three pairs of solid curves 
mark r = 0.3 (lower left) and 3 (upper right) for L49 = 0.1 (left), 1 (center), and 10 (right). For a 
given L, bursts that occur to the lower left of the r = 0.3 curve have afterglows that are relatively 



unextinguished by the circumburst cloud, and bursts that occur to the upper right of the r = 3 
curve have highly extinguished afterglows. Since circumburst extinction appears to be responsible 
for most of the dark bursts (Reichart & Yost 2001), we find in Reichart (2001b) that most of 
the dark bursts occur to the upper right of the r = 3 curve, which implies relatively large sizes 
and high densities for their circumburst clouds (§1). We show in Figures 1 and 2 that diffuse 
clouds of these sizes and densities are unlikely to be in equilibrium, or even near equilibrium, for 
typical temperatures and magnetic field strengths. Such clouds would be unstable to gravitational 
collapse, resulting in the collapse and fragmentation of the cloud until star formation re-establishes 
pressure equilibrium within the fragments, and the fragments are bound by self-gravity: This is 
precisely the structure of molecular clouds (e.g., Solomon et al. 1987). Furthermore, molecular 
clouds regularly have sizes R ;> 10 pc and masses M ;> 3 x 10 5 Mq, and are active regions 
of star formation (see §5.4). Consequently, dark bursts probably occur in collapsing diffuse 
clouds/forming molecular clouds, and/or molecular clouds that re-established equilibrium some 
time ago. We confirm that molecular clouds have sufficiently high column densities to make a 
large fraction of bursts dark at optical wavelengths in §3. 

Solomon et al. (1987) find the mass scale of the molecular cloud fragments to be on the order 
of a few solar masses, which supports the idea that they are supported against further collapse 
and fragmentation by star formation. Since the lifetimes of massive stars are comparable to the 
free fall time of the progenitor cloud for these cloud masses, the first generation of massive star 
formation and supernovae should occur on this short timescale: ~ 10 7 — 10 8 yr. Consequently, 
dark bursts are probably a byproduct of this burst of star formation if the molecular cloud formed 
recently, and/or the result of lingering or latter generation star formation if the molecular cloud 
formed some time ago. 



3. Consistency Check: The Column Density and Optical Depth Distributions of 

Bursts in Galactic-Like Molecular Clouds 

We now consider the column density and optical depth properties of bursts in Galactic-like 
molecular clouds.^ First, we compute the mean column densities of the 273 molecular clouds in 
the sample of Solomon et al. (1987), which was constructed using data from the Massachusetts- 
Stony Brook CO Galactic Plane Survey. The mean column density of a cloud is given by 
Nh = M/(fj,Hm}{A), where M is the virial mass (Solomon et al. 1987 show that the clouds are 
in or near virial equilibrium), A = 11. 56(D tan ^/aiOb) 2 is the effective cross section of the cloud 



We must entertain the possibility that Galactic molecular clouds are not representative of molecular clouds 
elsewhere in the universe, in which case the strength of our conclusions is diminished. However, since molecular clouds 
are formed from diffuse clouds that become too small or dense to support their mass, it seems that environmental 
differences between galaxies more likely affect the numbers of molecular clouds formed than the bulk physical 
properties of these clouds. This is supported by the uniformity of molecular cloud mean column densities across 
of the varied environments of our own galaxy (see below) . 



(Solomon et al. 1987), D is the distance to the cloud, and o\ and Ob are the angular sizes of the 
cloud. We plot the distribution (dotted histogram) of mean column densities in Figure 3, and 
confirm the finding of Solomon et al. (1987) that Ph^h is narrowly peaked around ~ 170 Mq 
pc~ 2 , independent of cloud mass and size. 

Next, we correct this distribution for a number of effects. First, since bursts more likely 
trace cloud mass than cloud number, we weight these mean column densities by cloud mass, and 
renormalize the distribution (dashed histogram). Since molecular cloud column densities appear 
to be fairly independent of cloud mass, the distribution is not significantly changed. Also, we 
weight the M < 7 x 10 4 M clouds by an additional factor of (M/7 x 10 4 M )~ 3 / 2 to correct 
for an undercounting of low mass clouds on the far side of the Galaxy (Solomon et al. 1987). 
However, this affects the mass- weighted mean column density distribution negligibly. 

Finally, we correct for geometrical effects: (1) bursts occur in their circumburst clouds, not 
behind them; and (2) molecular clouds are centrally condensed. First, we adopt the same density 
distribution for the clouds that Solomon et al. (1987) adopted to compute the virial masses: 
p oc r _1 within the effective radius of the cloud, and p = beyond this radius. This density 
distribution results in a surface brightness profile that is similar to what is observed (Solomon 
et al. 1987). Next, we trace this density distribution with 10 4 randomly-placed bursts for each 
cloud, and compute for each burst the distance and column density to the earth-facing side of the 
cloud. The effect of placing bursts in the clouds, as opposed to behind the clouds, is to lower the 
mean column densities on average by a factor of two, and sometimes (when bursts occur near 
the earth-facing side of the cloud) considerably more. However, this is partially offset by the 
central condensation of the clouds, which places more bursts in above average density (and column 
density) environments. We plot the final distribution - the expected column density distribution 
for bursts in Galactic- like molecular clouds - also in Figure 3 (solid histogram). 

Lastly, we confirm that Galactic-like molecular clouds, the column densities of which appear 
to be fairly universal, at least within the Galaxy, have sufficiently high column densities to make a 
large fraction of bursts dark at optical wavelengths. To this end, we plot in Figure 4 the expected 
distribution (1 and 2 a dashed contours) of column densities and distances to the earth-facing side 
of the cloud for bursts in Galactic-like molecular clouds. Also in Figure 4, we plot the curves of 
constant post-sublimation/fragmentation optical depth from Figures 1 and 2, and lines of constant 
column density. It appears that for low values of the isotropic-equivalent peak luminosity of the 
optical flash (L49 <J 0.1), Galactic- like molecular clouds can make a large fraction of bursts dark 
at optical wavelengths. However, even for very low values of L, a reasonable fraction of bursts 
would be left relatively unextinguished at optical wavelengths. Since comparable fractions of dark 
bursts and bursts with detected optical afterglows are observed (about 2/3 and 1/3, respectively; 
e.g., Fynbo et al. 2001; Lazzati, Covino & Ghisellini 2001), this suggests that the bursts with 
detected optical afterglows might also occur in molecular clouds. We confirm this by modeling 
and constraining the distribution of column densities, measured from absorption of the X-ray 
afterglow, of the bursts with detected optical afterglows in §4. We return to the issue of the 



isotropic-equivalent peak luminosity of the optical flash in §5.1. 



4. Evidence for a Molecular Cloud Origin for Bursts with Detected Optical 
Afterglows: The Column Density Distribution 

In this section, we model and constrain the column density distribution of the bursts 
with detected optical afterglows, and compare this distribution to the expected distribution 
for bursts in Galactic- like molecular clouds (§3). To this end, we have reviewed the literature, 
and list in Table 1 observer-frame column densities and uncertainties that have been measured 
from the X-ray afterglows of 15 bursts, 11 of which have detected optical afterglows (this 
information is not yet available for the vast majority of the dark bursts). Before modeling 
these data, we convert each measurement Njj,i and upper and lower 1-a uncertainty cr Ut i and 
an into probability distributions, first in the observer frame, and then in the source frame by 
correcting the observer-frame probability distribution for the Galactic column density along the 
line of sight, and the redshift of the burst if known. If a u % = an (or aij = Nh,i), we take the 
observer- frame probability distribution to be a Gaussian of mean Nn,i and standard deviation a u ,i'- 
Pi(NH,obs) = G(NH t0 b s ,NH t i,a u j). Otherwise, we take the observer-frame probability distribution 
to be a skewed Gaussian: Pi(Nu t obs) = G(Nfj obs , Nf Ii , af)dNfj obs /dNH,obs, where s is the skewness 
parameter, and s and <Tj are given by solving 

of = (N H ,i + v u ,i) s - N^ 

= N^ t -(N H ,-any, (4) 

i.e., we take the observer-frame probability distribution to be a Gaussian in Nfj obs , which reduces 
to a Gaussian in N}j,obs when a Ut i = ai t i (s = 1). The source- frame probability distribution 
Pi(Njj) is then given by substituting NH t obs = Nh,mw + Nh(1 + z)~ 2 ' 6 , where Nh,mw is the 
Galactic column density along the line of sight (interpolated^ from the maps of Dickey & Lockman 
1990), Nh is the source- frame column density, and (1 + z)~ 2 ^ scales Nh to the observer frame 
(e.g., Morrison & McCammon 1983). If multiple column density measurements are available for 
a burst, we average their source- frame probability distributions if the measurements were made 
from the same data, and we take the product of these distributions if the measurements were 
made from independent data (GRB 000926). We plot the source-frame probability distributions, 
peak-normalized to one, for the bursts with detected optical afterglows in the top panel of Figure 
5. The dotted curves mark bursts of unknown redshift, in which case we fix z = 0: These 
distributions can be thought of as fuzzy lower limits (see below). 

We model the column density distribution of the bursts with detected optical afterglows with 



We use the uh FTOOL at http://heasarc.gsfc.nasa.gov/cgi-bin/Tools/w3nh/w3nh.pl 



a broken power law: 

f (N H < 0) 

n(\ogN H ) = uho^E {^) b (0<N H <a) ; (5) 

IeTo^(^) C (N H >a) 

where a is the break column density, b is the low column density power-law index, and c is the high 
column density power-law index. We use a three-parameter function for flexibility, and a simple 
function (as opposed to, e.g., the more elegant skewed Gaussian) so the inner integral of Equation 
(P) can be evaluated analytically (otherwise, the error analysis would be computationally taxing). 
We note that Equation (||) is log-normalized [/f^, n(log Nn)dlog Nh = 1] to facilitate comparison 
of the fitted distribution to the solid histogram of Figure 3 (see below). 

We fit the model n(logiV#) to the data Pi(Nn) using Bayesian inference (e.g., Reichart 
2001a). The likelihood function is given by 

N 
C = ]j£i, (6) 

where N = 11 is the number of bursts with detected optical afterglows in our sample, and 

^{N^niN^dNu ( Zi known) 

(7) 
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(e.g., Reichart & Yost 2001), or equivalently 

IZoPi( N H)n(logN H )dlogN H [ Zi known) 
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since n(logiV^) = for N H < 0. We find that log(a/cm- 2 ) = 22.161^, arctan6 = 0.92±g;|| 
(6 w 1.31±g;||), and arctanc = — 1.46±g;ff [c ^ -1.33 (1 a)]. We plot the likelihood-weighted 
average model distribution (solid curve), which is also log-normalized, and the \-a uncertainty in 
this distribution (dotted curves) in the bottom panel of Figure 5. We find that the column density 
distribution of the bursts with detected optical afterglows peaks around 10 22 cm -2 , and spans a 
factor of a few to an order of magnitude or so. These findings verify and improve upon the earlier 
finding of Galama & Wijers (2001) that these column densities span 10 22 — 10 23 cm -2 . 

Lastly, we compare the fitted column density distribution of the bursts with detected optical 
afterglows to the expected column density distribution for bursts in Galactic-like molecular clouds 
(solid histogram, replotted from Figure 3), and find these distributions to be consistent within 
the l-o" uncertainties. Furthermore, the fitted distribution is not consistent with the expectation 
for bursts that occur in diffuse clouds, the canonical column density of which is a few times 10 20 
cm -2 . Consequently, we find that all but perhaps a few bursts, dark and otherwise, probably 
occur in molecular clouds. 



5. Discussion 

5.1. Implications for the Isotropic-Equivalent Peak Luminosity of the Optical Flash 

Now that we have shown that all but perhaps a few bursts probably occur in molecular 
clouds (§4), we return to the issue of the isotropic-equivalent peak luminosity of the optical flash. 
In §3, we showed that Galactic-like molecular clouds can make a large fraction of bursts dark 
at optical wavelengths, but only if a large fraction of bursts have L49 ^ 0.1. Furthermore, even 
if the value of L were very low, a reasonable fraction of these bursts would be sufficiently near 
the earth-facing side of their circumburst clouds, and/or in sufficiently small clouds, to be left 
relatively unextinguished at optical wavelengths (Figure 4). Since comparable fractions of dark 
bursts and bursts with detected optical afterglows are observed, this suggests that most bursts, 
dark or otherwise, have L49 <^ 0.1, and that bursts with L49 p> 1 should be relatively rare. Indeed, 
only one such burst, GRB 990123 (Akerlof et al. 1999) with L 49 « 1 (Waxman & Draine 2000), 
has been identified to date. Furthermore, of 11 other bursts that were observed on this timescale 
by ROTSE I, all have flux upper limits that are fainter than the peak flux of the optical flash of 
GRB 990123, often by 4 - 5 magnitudes (Akerlof et al. 2000; Kehoe et al. 2001). Similar results 
have been found with LOTIS (e.g., Williams et al. 2000). 

These findings are consistent with the simple model introduced by Reichart & Yost (2001), 
in which the collimation angle of a burst's ejecta, which strongly affects the value of L, and the 
column density to the earth-facing side of the circumburst cloud, determine whether a burst is 
dark or not. 



5.2. Implications for the Density of the Circumburst Medium 

Figure 4 suggests that bursts typically occur in regions of density n# ~ 30 — 10 3 cm~ 3 . 
However, to date, the modeling of X-ray through radio afterglow data has suggested significantly 
lower densities, ranging from n~5x 10~ 4 cm~ 3 (GRB 990123; Panaitescu & Kumar 2001) to 
n ~ 30 cm~ 3 (GRB 000926; Harrison et al. 2001). A number of reasonable solutions to this 
problem come to mind. First, the densities in Figure 4 are average densities to the earth- facing 
side of the circumburst cloud, but molecular clouds consist of fragments bound by self-gravity 
(§2). Solomon et al. (1987) estimate that there are ~ 10 fragments along a typical line of sight 
through a typical molecular cloud, that the density of these fragments is on the order of 10 3 cm~ 3 , 
and that the size of these fragments is on the order of a few tenths of a parsec, which is very 
roughly the distance a burst's ejecta travel before the burst occurs. Consequently, low densities 
might be expected if the ejecta clear the circumburst fragment before the burst occurs, or at least 



before the observed afterglow occurs.]] Alternatively, the circumburst fragment might be dispersed 
by prior massive star formation within the fragment (e.g., Scalo & Wheeler 2001). 

Another reasonable solution to this problem is that the practice of modeling afterglow data is 
not yet, or is just now becoming, sufficiently sophisticated to correctly extract the density of the 
circumburst medium. Case in point is that only very recently has the inverse-Compton component 
of the afterglow been included in modeling efforts: This component is essential to correctly 
interpret afterglow data if n ;> 1 cm~ 3 (e.g., Sari & Esin 2001). Including this component in their 
model, Harrison et al. (2001) found the first observational evidence for this component in the 
afterglow of GRB 000926: The strength of this component implied a density for the circumburst 
medium of n ~ 30 cm -3 , the highest value yet reported. However, more interesting is the fact that 
the source-frame column density of this burst is the most stringently low of all of the bursts (of 
known redshift) in the top panel of Figure 5. This suggests that even higher densities might be 
found for the circumburst environments of the other bursts with detected optical afterglows once 
properly modeled, or once this region of the parameter space is more fully explored. 

Lastly, we point out that these modeling efforts are naturally biased against the identification 
of high density circumburst environments. The self-absorption frequency increases with increasing 
density of the circumburst medium (e.g., Sari, Piran &; Narayan 1998; Sari & Esin 2001): Once 
n ^ 10 — 10 2 cm~ 3 (for typical values of other parameters; Reichart & Yost 2001), the radio 
afterglow becomes increasingly difficult to detect. Also, if a few parsecs of dust remain along the 
line of sight, once n ;> 10 2 cm -3 , the optical afterglow becomes increasing difficult to detect. Since 
modelers tend to avoid radio dark and/or optically dark bursts, they naturally bias themselves 
against the identification of circumburst environments of these densities. 

We address this problem in a simple way in Reichart & Yost (2001): We show that a range of 
densities for the circumburst medium, with other parameters fixed to canonical values, can explain 
snapshot optical vs. X-ray and optical vs. radio brightness distributions of the afterglows of all 
bursts, radio dark and optically dark bursts included. We find that the density of the circumburst 
medium probably spans many orders of magnitude, from densities possibly as low as densities that 
are typical of the Galactic disk to densities probably as high as densities that are typical of dense 
clouds. 



7 If the ejecta clear the circumburst fragment after the burst occurs, one might be able to measure a change in 
the X-ray column density. Although this would be only a f - 20% effect for most bursts (assuming typically 5 - f 
fragments between the burst and the earth- facing side of the cloud), a few bursts should occur sufficiently close to 
the earth-facing side of the cloud to make this a 50 - 100% effect. 
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5.3. Implications for Circumburst Extinction Studies 

Figure 4 also suggests that if most bursts have L49 <J 1, which appears to be the case (§5.1), 
a reasonable fraction of bursts should have afterglows that are extinguished at intermediate levels: 
0.3 <J r <J 3, where r ~ source-frame Ab (Reichart 2001b). Such bursts would be the missing 
link between dark bursts and bursts with relatively unextinguished afterglows. Furthermore, by 
extracting extinction curves from the afterglow data of such bursts (Reichart 2001a), bursts could 
serve as probes of dust under the extreme conditions of sublimation and fragmentation, and as 
probes of dust as a function of redshift, to very high redshifts (Lamb & Reichart 2000). 



5.4. Implications for the Nature of Star Formation in the Universe 

Mooney & Solomon (1988) find that the far- infrared luminosity-to-mass ratio for isolated and 
weakly interacting spiral galaxies is consistent with the average for Galactic molecular clouds, 
which suggests that most, if not all, of the star formation in such galaxies occurs in molecular 
clouds. Our finding that most, if not all, bursts occur in molecular clouds would strengthen this 
idea if most of the star formation in the universe occurs in such galaxies (see Lamb & Reichart 
2000 for a review of evidence linking bursts to star formation and massive star death) . 

However, Ramirez-Ruiz, Trentham & Blain (2001) argue that about 3/4 of the star formation 
in the universe occurs in ultraluminous infrared/submillimeter-bright galaxies. Often, this emission 
is concentrated within hundreds of parsecs of the center of the galaxy. Since the column densities 
to stars in these nuclear regions range from Njj ~ 10 23 — 10 24 cm~ 2 , and the sizes of these regions 
are one to two orders of magnitude too large for dust sublimation and fragmentation to matter, 
bursts that occur in these regions would necessarily be dark at optical wavelengths. Since about 
2/3 of bursts are dark at optical wavelengths, Ramirez-Ruiz, Trentham & Blain (2001) argue that 
most of the dark bursts probably occur in such regions. 

We now test this idea. In Figure 6, we plot the source- frame column density probability 
distributions of the four dark bursts in Table 1 (§4). Redshifts are not known for three of these 
bursts, so we adopt z = 1 (top panel), 2 (middle panel), and 3 (bottom panel) for these three 
bursts. Given the redshifts that have been measured for bursts to date, the probabilities that all 
three of these bursts have z > 1, 2, and 3 are 0.27, 0.013, and 0.0016, respectively. Even in the 
very unlikely event that all three of these bursts have z ~ 3, these four probability distributions 
are not consistent with the expected column density distribution for bursts that occur in the 
nuclear regions of ultraluminous infrared/submillimeter-bright galaxies, but are consistent with 
the expected column density distribution for bursts that occur in molecular clouds (§§3 and 4). 

However, these ideas can be reconciled if the stars, and in particular the massive stars, form 
predominately outside of the nuclear region. Ultraluminous infrared/submillimeter-bright galaxies 
are thought to be strongly interacting spiral galaxies, where the interaction causes a starburst 
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(e.g., Solomon & Sage 1988; Solomon et al. 1997). The light from the starburst is absorbed by 
dust in the star-forming regions, and is re-emitted in the far infrared. We now consider how the 
interaction causes the starburst, and where the stars are when this occurs. We consider first 
the early stages of the interaction. Solomon et al. (1997) argue that the starburst cannot be 
caused by the collisions of molecular clouds, or by the creation of new molecular clouds from the 
collisions of diffuse clouds, since the interaction velocity exceeds the escape velocity of even the 
most massive clouds by more than an order of magnitude. Furthermore, they argue that molecular 
cloud collisions are improbable, since their mean free path exceeds typical disk thicknesses by 
about an order of magnitude. However, diffuse cloud collisions are very probable, the result of 
which should be a hot, ionized gas that surrounds the molecular clouds in that region of the 
galaxy (Jog & Solomon 1992). This high-pressure gas should cause radiative shock compression 
of the outer layers of these molecular clouds, resulting in a starburst for which Jog & Solomon 
(1992) estimate a far-infrared luminosity of a few times 10 11 Lq. Starbursts of this nature would 
span large regions of the galaxy, and indeed, starbursts in colliding galaxies are observed to span 
several kiloparsecs (e.g., Wright et al. 1984; Joseph & Wright 1985; Rieke et al. 1985; Wright et 
al. 1988). Such starbursts are consistent with Figure 6. 

We now consider the latter stages of the interaction. Tidal interaction is thought to create a 
dense nuclear region: The perturbation of cloud orbits increases the diffuse cloud collision rate, 
which results in a dissipation of kinetic energy that causes most of the gas in the disk to fall to 
the center of the galaxy (e.g., Noguchi & Ishibashi 1986; Hernquist 1989). However, dense gas 
alone does not a starburst make. Jog & Das (1992) propose that nuclear starbursts are caused by 
molecular clouds that tumble into the nuclear region: Again, high-pressure gas causes radiative 
shock compression of the outer layers of the molecular clouds, resulting in a starburst for which 
Jog & Das (1992) estimate a far- infrared luminosity that ranges from ;> 10 9 Lq for distant tidal 
interactions to <^ 10 12 Lq for very perturbed mergers. Even nuclear starbursts might be consistent 
with Figure 6 if the molecular clouds are shock compressed ~ 10 7 — 10 8 yr before they fall into 
the very high column density environment of the inner nuclear region: This would give them 
sufficient time to form and expend their massive stars in low and/or intermediate column density 
environments. 

Alternatively, the fraction of star formation in the universe that occurs in ultraluminous 
infrared/submillimeter-bright galaxies might be smaller than ~ 3/4. This would imply that (1) the 
IMF is high-massed biased in SCUBA galaxies, (2) SCUBA sources are powered predominately by 
AGN, and/or (3) the star-formation rates in UV-bright galaxies are underestimated (Ramirez- Ruiz, 
Trentham & Blain 2001). 

However, the idea that bursts should be associated with ultraluminous infrared/submillimeter- 
bright galaxies is supported by the recent associations of GRB 980703 with an ultraluminous 
infrared galaxy (Berger, Kulkarni & Frail 2001), and GRB 010222 with a submillimeter-bright 
galaxy (Frail et al. 2001). However, since both GRB 980703 and GRB 010222 have well-detected 
optical afterglows, they could not have occurred in dense nuclear regions, and hence are in 
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agreement with Figure 6 and the above possible explanations. 

In any case, our finding that most, if not all, bursts occur in molecular clouds strengthens the 
idea that most of the star formation in the universe occurs in molecular clouds, be they in isolated 
or interacting galaxies. 



6. Conclusions 

We find that the circumburst clouds of dark bursts are too large and massive to be diffuse, and 
consequently are probably molecular, making them active regions of star formation. Furthermore, 
we show that Galactic molecular clouds, the column densities of which appear to be fairly 
universal, at least within the Galaxy, have sufficiently high column densities to make a large 
fraction of bursts dark at optical wavelengths. However, since the optical flash sublimates dust 
and the burst and afterglow fragment dust, we find that most bursts in Galactic- like molecular 
clouds would have to have L49 ^0.1: This seems reasonable, given the relatively deep imaging 
of ROTSE I and LOTIS. Furthermore, we find that even for very low values of L, a reasonable 
fraction of bursts in Galactic-like molecular clouds would be sufficiently near the earth-facing 
side of the cloud, and/or in sufficiently small clouds, to be left relatively unextinguished at 
optical wavelengths. Consequently, we model and constrain the distribution of column densities, 
measured from absorption of the X-ray afterglow, of the bursts with detected optical afterglows: 
We find that these bursts are also consistent with a molecular cloud origin, and are not consistent 
with a diffuse cloud origin. Consequently, we find that all but perhaps a few bursts, dark or 
otherwise, probably occur in molecular clouds: This seems reasonable, given that it is thought 
that most of the star formation in the universe occurs in molecular clouds, be they in isolated or 
interacting galaxies. Finally, we show that the limited information that is available on the column 
densities of the dark bursts is not consistent with the idea that the dark bursts occur in the 
nuclear regions of ultraluminous infrared/submillimeter-bright galaxies. However, this does not 
necessarily imply that most of the star formation in the universe does not occur in ultraluminous 
infrared/submillimeter-bright galaxies. 
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Table 1. Column Density Data 



GRB 


z 


Dark 


JVff,Mlf a ' 


NH,obs a 


ReF 


970228 
970402 


0.695 


N 
Y 


1.60 
2.04 


o c+3.3 
°-°-2.3 

65+ 30 
U - DO -0.65 


1 

2 
2 


970508 


0.835 


N 


0.516 


1.8+fl 


2 


970828 
971214 


0.958 
3.418 


Y 

N 


0.356 
0.162 


q c+2.2 

4-i±?:S 

4.3 + |^ 

1+2.3 
L -l 


2 
3 

2 
4 


980329 
980519 




N 
N 


0.941 

1.74 


10 ±4 
10.71U 

5.6^ 2 


5 
2 
2 


990703 


0.966 


N 


0.580 


6.6±i;l d 


6 


990123 
990510 


1.600 
1.619 


N 
N 


0.212 
0.937 


i rr+0.7 
l -°-0.6 

2.1 ±0.6 
1.7 ±0.7 
1.5 ±0.6 


7 
8 
8 
8 


991014 




Y 


2.45 


< 8.6 


9 


991216 


1.020 


N 


2.04 


3.5 ±1.5 


10 


000214 




Y 


0.577 


n n? +0.75 


11 


000926 


2.037 


N 


0.266 


0.5 ±0.3 
0.48 ±0.3 
0.30 ± 0.25 c 


12 
13 
13 


010222 


1.477 


N 


0.162 


1.5 ±0.3 


14 



a 10 21 cm' 2 . 

Interpolated from the maps of Dickey & Lockman 
1990. 

c l. Frontera et al. 1998; 2. Owens et al. 1998; 3. 
Yoshida et al. 1999; 4. Dal Fiume et al. 2000; 5. in 't 
Zand et al. 1998; 6. Vreeswijk et al. 1999; 7. Yonetoku 
et al. 2000; 8. Kuulkers et al. 2000; 9. in 't Zand et al. 
2000; 10. Piro et al. 2000; 11. Antonelli et al. 2000; 12. 
Piro et al. 2001; 13. Harrison et al. 2001; 14. in 't Zand 
et al. 2001. 

Determined from information in Vreeswijk et al. 
1999. 

c Measurement from independent data. 
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Fig. 1. — Post-sublimation/fragmentation optical depth r to a burst that is embedded a distance r 
within a cloud of constant hydrogen density hh, along the line of sight (solid curves; from Reichart 
2001b), and diffuse cloud stability criteria (dashed curves). The three pairs of solid curves mark 
r = 0.3 (lower left) and 3 (upper right) for L49 = 0.1 (left), 1 (center), and 10 (right). For a 
given L, bursts that occur to the lower left of the r = 0.3 curve have afterglows that are relatively 
unextinguished by the circumburst cloud, and bursts that occur to the upper right of the r = 3 
curve have highly extinguished afterglows. The dashed curves mark the maximum radius for which 
a diffuse cloud is stable against gravitational collapse, for cloud temperatures T = 10 and 100 K 
and cloud magnetic field strength B = 3 /iG (§2). 

Fig. 2. — Same as Figure 1, except for T = 30 K and B = 1 and 10 fiG (§2). 

Fig. 3. — Expected column density distribution n(\ogNjj) (log-normalized to one) for bursts that 
occur in Galactic-like molecular clouds (solid histogram). The dotted histogram marks the mean 
column density distribution of the Galactic molecular clouds of Solomon et al. (1987), which we 
compute from their mass and size measurements. The dashed histogram marks the mass-weighted 
mean column density distribution, which we also correct for an undercounting of low mass clouds 
on the far side of the Galaxy (a negligible correction). The solid histogram marks this distribution 
after correcting for the fact that bursts occur in their circumburst clouds, not behind them, and 
for the fact that molecular clouds are centrally condensed (§3). 

Fig. 4. — Expected distribution of column densities and distances to the earth-facing side of the 
cloud for bursts that occur in Galactic- like molecular clouds (dashed contours), and the curves of 
constant post-sublimation/fragmentation optical depth from Figures 1 and 2 (solid curves). The 
thick dashed contour contains 68.3% of the distribution, and the thin dashed contour contains 
95.4% of the distribution. The dotted lines mark constant hydrogen column densities (§3). 

Fig. 5. — Top panel: Source-frame column density probability distributions Pi{Nf{) (peak- 
normalized to one) for the 11 bursts with detected optical afterglows in Table 1. The dotted curves 
mark bursts of unknown redshift, for which we fix z = 0: These distributions can be thought of 
as fuzzy lower limits. Bottom panel: Likelihood- weighted average model distribution n(logA^jy) 
(solid curve), and 1-a uncertainty in this distribution (dotted curves). The solid histogram is the 
expected column density distribution for bursts that occur in Galactic-like molecular clouds from 
Figure 3. These distributions are consistent within the 1-a uncertainties (§4). 

Fig. 6. — Source-frame column density probability distributions Pi{Nh) for the four dark bursts 
in Table 1. The dotted curves mark bursts of unknown redshift, for which we adopt z = 1 (top 
panel), 2 (middle panel), and 3 (bottom panel; §5.4). 
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